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Abstract 

The recently-developed non-equilibrium extension of the self-consistent generalized Langevin 
equation theory of irreversible relaxation [Phys. Rev. E (2010) 82, 061503; ibid. 061504] is applied 
to the description of the irreversible process of equilibration and aging of a glass-forming soft-sphere 
liquid that follows a sudden temperature quench, within the constraint that the local mean particle 
density remains uniform and constant. For these particular conditions, this theory describes the 
non-equilibrium evolution of the static structure factor S(k; t) and of the dynamic properties, such 
as the self-intermediate scattering function Fs(k, r; t), where r is the correlation delay time and t is 
the evolution or waiting time after the quench. Specific predictions are presented, for the deepest 
quench (to zero temperature). The predicted evolution of the o-relaxation time r a (t) as a function 
of t allows us to define the equilibration time t eq ((p), as the time after which r a (t) has attained 
its equilibrium value Ta q ((f>). It is predicted that both, t eq {4>) and Ta 9 ((ft), diverge as <p — > <p^ a \ 
where (p^ is the hard-sphere dynamic-arrest volume fraction (p^ (~ 0.582), thus suggesting that 
the measurement of equilibrium properties at and above (jr- a > is experimentally impossible. The 
theory also predicts that for fixed finite waiting times t, the plot of T a (t;4>) as a function of <j) 
exhibits two regimes, corresponding to samples that have fully equilibrated within this waiting 
time ((f) < <f)( c '(t)), and to samples for which equilibration is not yet complete (4> > (f>^ c \t)). The 
crossover volume fraction (t) increases with t but saturates to the value cf>^ . 

PACS numbers: 23.23,+x, 56.65.Dy 
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I. INTRODUCTION. 

Classical and statistical thermodynamics deal with the equilibrium states of matter jl, 2]. 
Driving the system from one equilibrium state to another, however, involves the passage of 
the system through a sequence of instantaneous states that do not satisfy the conditions 
for thermodynamic equilibrium, and hence, constitute a non-equilibrium process 3|, |4|. The 
description of these processes fall outside the realm of classical and statistical thermodynam- 
ics, unless the sequence of non-equilibrium states do not depart appreciably from a sequence 
of equilibrium states. Such idealized process can be thought of as an infinite sequence of 
infinitesimally small changes in the driving control parameter, after each of which the system 
is given sufficient time to equilibrate. This so-called quasistatic process is an excellent rep- 
resentation of real process when the equilibration times of the system are sufficiently short. 
However, when the equilibration kinetics is very slow, virtually any change will involve in- 
trinsically non-equilibrium states whose fundamental understanding must unavoidably be 
done from the perspective of a non-equilibrium theory ja]. 

These concepts become particularly relevant for the description of the slow dynamics 
of metastable glass-forming liquids in the vicinity of the glass transition 

aa. 

It is well 

known that the decay time of the slowest relaxation processes (the so-called a-relaxation 
time r a ) increases without bound as the temperature T is lowered below the glass transition 
temperature T^ 9 \ It is then natural to think that the equilibration time of the system 
must also increase accordingly. To be more precise, let us imagine that a glass- forming 
liquid, initially at an arbitrary temperature is suddenly cooled at time t = to a final 
temperature T, after which it is allowed to evolve spontaneously toward its thermodynamic 
equilibrium state. Imagine that we then monitor its a- relaxation time r a (t) as a function 
of the evolution or "waiting" time t elapsed after the quench. We say that the system has 
equilibrated when r a (t) reaches the plateau that defines its final equilibrium value t^ q (T), 
which must only depend on the final temperature T. The beginning of this plateau occurs 
at a certain value of the waiting time t, that we refer to as the equilibration time t eq (T); this 
equilibration time must also depend on the final temperature T. 

There are strong indications, from recent computer simulation experiments , that in 
the metastable regime these two characteristic times, t^{T) and t eq (T), are related to each 
other as t eq (T) oc [T^ q (T)] v , with an exponent r] > 1 (more specifically, rj w 1.5 [si, This 
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implies that in order to measure the actual equilibrium value t^(T) we have to wait, before 
starting the measurement of r^ q (T), for an equilibration time t eq (T) that will increase faster 
than t^ q (T) itself. This poses an obvious practical problem for the measurement of r^ q (T) 
when the temperature T approaches the glass transition temperature T^ 9 \ since sooner or 
later we shall be unable to wait this required equilibration time. This situation then implies 
that it is impossible to discard a scenario in which the equilibrium a-relaxation time T^ q (T) 
diverges at a singular temperature T^ a \ since the equilibration time t eg (T) needed to observe 
this divergence will also diverge at that temperature, i.e., it will be impossible to equilibrate 
the system at a final temperature near or below within experimental waiting times. Of 
course, a measurement carried out at a finite t, will always report a result for r a , but this 
result will correspond to r a (t), the non-equilibrium value of the a-relaxation time registered 
at that waiting time t. Thus, the analysis of these experimental measurements cannot 
be based on the postulate that the system has reached equilibrium; instead, one needs to 
interpret these experiments in the framework of a quantitative theory of slowly-relaxing 
non-equilibrium processes. 

Until recently, however, no quantitative, first-principles theory had been developed and 
applied to describe the slow non- equilibrium relaxation of structural glass-forming atomic 
or colloidal liquids. About a decade ago Latz 10[ attempted to extend the conventional 
mode coupling theory (MCT) of the ideal glass transition llMl4J|. to describe the aging of 
suddenly quenched glass forming liquids. A major aspect of his work involved the gener- 
alization to non-equilibrium conditions of the conventional equilibrium projection operator 
approach 15| to derive the corresponding memory function equations in which the mode 
coupling approximations could be introduced. Similarly, De Gregorio et al. [if]] discussed 
time-translational invariance and the fluctuation-dissipation theorem in the context of the 
description of slow dynamics in system out of equilibrium but close to dynamical arrest. 
They also proposed extensions of approximations long known within MCT. Unfortunately, 
in neither of these theoretical efforts, quantitative predictions were presented that could be 
contrasted with experimental or simulated results in specific model systems of structural 
glass-formers. 

In an independent but similarly-aimed effort, on the other hand, the self-consistent gen- 
eralized Langevin equation (SCGLE) theory of colloid dynamics 17H20| and of dynamic 
arrest (21 25] has recently been extended to describe the (non-equilibrium) spatially non- 



3 



uniform and temporally non-stationary evolution of glass-forming colloidal liquids. Such an 



extension was introduced and described in detail in Ref. [26j, and will be referred to as the 
non- equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory. As one 
can imagine, the number and variety of the phenomena that could be studied with this new 
theory may be enormous, and to start its systematic application we must focus on simple 
classes of physically relevant conditions. Thus, as a first simple illustrative application, this 



theory was applied in Ref. 



271 ] to a model colloidal liquid with hard-sphere plus short-ranged 



attractive interactions, suddenly quenched to an attractive glass state. 

The aim of the present work is to start a systematic exploration of the scenario predicted 
by this theory when applied to the simplest irreversible processes, in the simplest and best- 
defined model system. In the present case we refer to the irreversible isochoric evolution of a 
glass-forming liquid of particles interacting through purely repulsive soft-sphere interactions, 
initially at a fluid-like state, whose temperature is suddenly quenched to a final value T^> = 
0, at which the expected equilibrium state is that of a hard-sphere liquid at volume fraction 
<p. Such process mimics the spontaneous search for the equilibrium state of this hard- 
sphere liquid, driven to non-equilibrium conditions by some perturbation (shear, for example 



29j ) which ceases at a time t = 0. One possibility is that the system will recover its 



equilibrium state within an equilibration time t eq (<p) that depends on the fixed volume 
fraction <p. The other possibility is that the system ages forever in the process of becoming 
a glass. The application of the NE-SCGLE theory to these irreversible processes results in 
a well-defined scenario of the spontaneous non-equilibrium response of the system, whose 
main features are explained and illustrated in this paper. 

In the following section we provide a brief summary of the non-equilibrium self-consistent 
generalized Langevin equation theory, appropriately written to describe the equilibration of 
a monocomponent glass-forming liquid constrained to remain spatially uniform. Section II III 
defines the specific model to which this theory will be applied, discusses the strategy of so- 
lution of the resulting equations, and illustrates the main features of the results. Section [TV] 
presents the scenario predicted by the NE-SCGLE theory for the first possibility mentioned 
above, namely, that the system is able to reach its thermodynamic equilibrium state. In this 
case we find that the equilibrium a-relaxation time t^ 9 (0), and the equilibration time t eq (4>) 
needed to reach it, will remain finite for volume fractions smaller than a critical value cf)^ a \ but 
that both characteristic times will diverge as <fi approaches this dynamic-arrest volume frac- 
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tion <j)^ « 0.582, and will remain infinite for > <p^ a \ Although it is intrinsically impossible 
to witness the actual predicted divergence, the theory makes distinct predictions regarding 
the transient non-equilibrium evolution occurring within experimentally-reasonable waiting 
times t. 

In Sect. |V]we analyze the complementary regime, > <p^ a \ in which the system, rather 
than reaching equilibrium within finite waiting times, is predicted to age forever. In this 
regime we find that the long-time asymptotic limit of S(k; t) will no longer be the expected 
equilibrium static structure factor S^ eq \k), but another, non-equilibrium but well-defined, 
static structure factor, that we denote as S^ a \k), and which depends on the protocol of 
the quench. Furthermore, contrary to the kinetics of the equilibration process, in which 
S(k;t) approaches S^ eg \k) in an exponential-like fashion, this time the decay of S(k;t) to 
its asymptotic value S^(k) follows a much slower power law. 

In section [Vj we put together the two regimes just described, in an integrated picture, 
which outlines the predicted scenario for the crossover from equilibration to aging. There 
we find that the discontinuous and singular behavior underlying the previous scenario is 
intrinsically unobservable, due to the finiteness of the experimental measurements, which 
constraints the observations to finite time windows. This practical but fundamental limita- 
tion converts the discontinuous dynamic arrest transition into a blurred crossover, strongly 
dependent on the protocol of the experiment and of the measurements. 

The main purpose of the present paper is to explain in sufficient detail the methodolog- 
ical aspects of the application of the theory, so as to serve as a reliable reference for the 
eventual application of this non-equilibrium theory to the same system but with different 
non-equilibrium processes (e.g., different quench protocols), or in general to different sys- 
tems and processes. Thus, we shall not report here the results of the systematic quantitative 
comparison of the scenario explained here with available specific simulations or experiments, 
which are being reported separately. Thus, the final section of the paper briefly refers to 
the main features of those comparisons, and discusses possible directions for further work. 

II. REVIEW OF THE NE-SCGLE THEORY. 

Let us mention that the referred non-equilibrium self-consistent generalized Langevin 
equation (NE-SCGLE) theory derives from a non-equilibrium extension of Onsager's theory 
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of thermal fluctuations [26(, and it consists of the time evolution equations for the mean 



value n(r, t) and for the covariance cx(r, r'; t) = 8n(r, t)Sn(r\ t) of the fluctuations 5n(r, t) = 
n(r,t) — n(r,t) of the local concentration profile n(r,t) of a colloidal liquid. These two 
equations are coupled, through a local mobility function b(r, t), with the two-time correlation 



function C(r, r'; t, t') = <5n(r, t)5n(r', t'). A set of well-defined approximations on the memory 
function of C(r, r'; t, t'), detailed in Ref. 26], results in the referred NE-SCGLE theory. 

As discussed in Ref. j^fj], for given interparticle interactions and applied external fields, 
the NE-SCGLE self-consistent theory is in principle able to describe the evolution of a 
strongly correlated liquid from an initial state with arbitrary mean and covariance n°(r) 
and cr°(r, r'), towards its equilibrium state characterized by the equilibrium local concen- 
tration profile n eq (r) and equilibrium covariance a eq (r, r'). These equations are in principle 
quite general, and contain well known theories as particular limits. For example, ignoring 
certain memory function effects, the evolution equation for the mean profile n(r, t) becomes 
the fundamental equation of dynamic density functional theory jsoj], whereas the "conven- 
tional" equilibrium SCGLE theory [23[ (analogous in most senses to MCT 11] ) is recovered 
when full equilibration is assumed and spatial heterogeneities are suppressed. The NE- 
SCGLE theory, however, provides a much more general theoretical framework, which in 
principle describes the spatially heterogeneous and temporally non- stationary evolution of a 
liquid toward its ordinary stable thermodynamic equilibrium state. This state, however, will 
become unreachable if well-defined dynamic arrest conditions arise along the equilibration 
pathway, in which case the system evolves towards a distinct and predictable dynamically 
arrested state through an evolution process that involves aging as an essential feature. 

To start the systematic application of this general theory to more specific phenomena we 
must focus on a simple class of physical conditions. Thus, let us consider the irreversible 
evolution of the structure and dynamics of a system constrained to suffer a programmed 
process of spatially homogeneous compression or expansion (and/or of cooling or heating). 
Under these conditions, rather than solving the time-evolution equation for n(r;t), we as- 
sume that the system is constrained to remain spatially uniform, n(r; t) = n(t), according to 
a prescribed time-dependence n{t) of the uniform bulk concentration and/or to a prescribed 
uniform time- dependent temperature T(t). Among the many possible programmed protocols 
(n(t), T(t)) that one could devise to drive or to prepare the system, in this paper we restrict 
ourselves to one of the simplest and most fundamental protocols, which corresponds to the 
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limit in which the system, initially at an equilibrium state determined by initial values of the 
control parameters, (n^\T^'), must adjust itself in response to a sudden and instantaneous 
change of these control parameters to new values (n^' ,T^'), according to the "program" 
n(t) = n^6(-t) + n^e(t) and T(t) = T^9(-t)+T^9(t), with 9{t) being Heavyside's step 
function. Furthermore, just like in the first illustrative example described in Ref. j^], here 
we shall also restrict ourselves to the description of an even simpler subclass of irreversible 
processes, namely, the isochoric cooling or heating of the system, in which its number density 
is constrained to remain constant, i.e., n(t) = n^' = rffl = n, while the temperature T(t) 
changes abruptly from its initial constant value T^' to a final constant value at t — 0. 
Under conditions of spatial uniformity, C(r, r'; t, t') can be written as 

C(\v-T'\,t'-t; t) = —— 3 J dkexp[-ik-(v-v')}F(k,T;t), (2.1) 

with t = {f — t) > 0, and where F(k, r; t) is the ^-evolving non-equilibrium intermediate 
scattering function (NE-ISF). Similarly, the covariance <r(r, r';t) can be written as 

?T f 

a(|r-r'|;i) = — J dkexpf-ik • (r - r')]^; t) (2.2) 

with S(k; t) = F(k, r = 0; t) being the time-evolving static structure factor. Under these 
conditions, the NE-SCGLE theory determines that the time-evolution equation for the co- 

27l |) may be written as an equation for S(k; t) which, for t > 0, 



variance (Eq. (2.11) of Ref. 
reads 

= -2k 2 D%(t)n^£^\k) [S(k; t) - l/n£«\k)] . (2.3) 

In this equation the function £V>(k) = £(k;n,T^) is the Fourier transform (FT) of the 
functional derivative £[\ r — r' \\n,T] = [5/3//[r; n\/8n(r')], evaluated at n(r) = n and T = 
As discussed in Refs. 26|, |27j, this thermodynamic object embodies the information, 



assumed known, of the chemical equation of state, i.e., of the functional dependence of the 
electrochemical potential /i[r;n] on the number density profile n(r). 

The solution of this equation, for arbitrary initial condition S(k; t = 0) = S^ l '(k), can be 
written as 

S(k; t) = S {i \k)e~ a{k)u{t) + [n^ik)}- 1 (l - e - a{k)u{t) ) , (2.4) 

with 

a{k) = 2k 2 D°n£ {f \k), (2.5) 
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and with 



u(t) = / b(t')dt'. (2.6) 
Jo 

In the equations above, the time-evolving mobility b(t) is defined as b(t) = DL(t)/D , 
with D° being the short-time self-diffusion coefficient and Dj,{t) the long-time self-diffusion 
coefficient at evolution time t. As explained in Refs. {26] and [27]], the equation 

/■oo 

b(t) = [l+ / drACi^t)]- 1 (2.7) 
Jo 

relates b(t) with the t-evolving, r-dependent friction coefficient A(*(t; t) given approximately 
by 



S(k;t) - 1 
S{k;t) 



(2.8) 



xF(k,r;t)F s (k,r;t). 

Thus, the presence of b(t) in Eq. ( 12. 6ft couples the formal solution for 5*(A;; t) in Eq. (12.4ft 
with the solution of the non-equilibrium version of the SCGLE equations for the collective 
and self NE-ISFs F(k,r;t) and Fs(k,z;t). These equations are written, in terms of the 
Laplace transforms (LT) F(k,z;t) and F s (k,r;t), as 

F(k, z;t) — k^pos-^(k;t) ' ( 2 - 9 ) 

Z 1+A(fc) AC*(z;t) 

and 

F 5 (fc, z- 1) = \ 2D0 , (2.10) 

Z + 1+A(fc) AC*(2;t) 

with X(k) being a phenomenological "interpolating function" [23], given by 

X(k) = 1/[1 + (A;/A; c ) 2 ], (2.11) 

with fc c = 1.305 x k max (t), where k max (t) is the position of the main peak of S(k;t) (in 
practice, however, k c w 1.305 (2tt/(t) [8]). The simultaneous solution of Ecs. (12. 3|) - (12. 10[) 
above, constitute the NE-SCGLE description of the spontaneous evolution of the structure 
and dynamics of an instantaneously and homogeneously quenched liquid. 

Of course, one important aspect of this analysis refers to the possibility that along the 
process the system happens to reach the condition of dynamic arrest. For the discussion of 
this important aspect it is useful to consider the long-r (or small z) asymptotic stationary 
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solutions of Eqs. (I2.9p - (l2.8p . the so-called non-ergodicity parameters, which are given by 



26] 



and 



W'V-gL S(k) " X(k;t)S(k;t) + k^(t) [ZAZ) 

/ g (fc;t) S ] teW ,r;t) = A(fc; ^g >7(t) > (2.13) 
where the t-dependent squared localization length 7(2) is the solution of 

1 (2-14) 



j(t) 67c 2 n^ Jo [X(k; t)S(k; t) + k 2 j(t)} [X(k; t) + k*<y{t)] ' 

Notice also that these equations are the non-equilibrium extension of the corresponding 
results of the equilibrium SCGLE theory (referred to as the "bifurcation equations" in 



the context of MCT 



111]), and their derivation from Eqs. f !2.8j) — (I2.10p follows the same 



arguments as in the equilibrium case [18] . The solution <y(t) of Eq. (I2.14p and the mobility 
b(t) constitute two complementary dynamic order parameters, in the sense that if ^(t) is 
finite (or b(t) = 0), then the system must be considered dynamically arrested at that waiting 
time t, whereas if j(t) is infinite, then the particles retain a finite mobility, b(t) > 0, and 
the instantaneous state of the system is ergodic or fluid-like. 

We recall that the first relevant application of Eq. (I2.14p is the determination of the 
equilibrium dynamic arrest diagram in control-parameter space (which, in the present case, 
is the density-temperature plane (n,T)). This diagram determines the region of fluid-like 
states, for which the solution ^ eq (n, T) (of Eq. f )2.14p . with S(k; t) = S eq (k; n, T)) is infinite. 
The complementary region contains the dynamically-arrested states, for which / y eq (n,T) is 
finite. The borderline between these two regions is the dynamic arrest transition line. Due 
to the complementarity of the dynamic order parameters j(t) and b(t), this curve is also 
the borderline between the region where the mobility b(t) will reach its equilibrium value, 
lim t _j. 0O b{t) = b eq (n, T) > 0, and the region of arrested states, where lim^oo b(t) = 0. Thus, 
since b eq (n, T) = D*(n, T) = Dl(u, T)/D°, where D 1,(71, T) is the equilibrium long-time self- 
diffusion coefficient at the point (n, T), this line is also the iso-diffusivity curve corresponding 
to D* = 0. 
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III. GENERAL FEATURES OF THE SOLUTION AND A SPECIFIC ILLUSTRA- 
TION. 



Let us now discuss some general features of the solution of the NE-SCGLE equations 
just presented. This discussion has a general character, but for the sake of clarity we shall 
illustrate the main concepts in the context of one specific application. Thus, consider a 
mono-component fluid of soft spheres of diameter a, whose particles interact through the 
truncated Lennard- Jones (TLJ) pair potential that vanishes for r > a, but which for r < a 
is given, in units of the thermal energy k B T = (3 , by 



(3u(r) 



r ) \r . 



The state space of this system is spanned by the volume fraction 
reduced temperature T* = k B T/e. 



(3.1) 

una 3 / 6 and the 



A. Thermodynamic framework: local curvature of the free energy surface. 

In order to apply Eqs. fl2.3l) - fl2.10l) to this model system, we first need to determine its 
thermodynamic property £^'(k). As indicated above, this is the Fourier transform of the 
functional derivative £[\ r— r' |; n, T] = [5/3/i[r; n]/5n(r')], which can also be written as £[\ r— 
r' \;n, T] = 8(r — r')/n — c(| r — r' |;n,T),with c(r; n, T) being the ordinary direct correlation 
function [2j. This is an intrinsically thermodynamic property, related with the equilibrium 
static structure factor S^ eq \k;n,T) by the Ornstein-Zernike (OZ) equation, which in Fourier 
space reads n£(k; n, T)S^ eq \k] n, T) = 1. The OZ equation is the basis for the construction of 
the approximate integral equations of the equilibrium statistical thermodynamics of liquids 
(2]]. In fact, we shall employ one such approximation to determine n£(k;n,T) for our soft- 
sphere system. This approximation, explained in detail in the appendix of Ref. [3l| and 
denoted as PY/VW, is based on the Percus-Yevick approximation [32[ within the Verlet- 
Weis correction |33| for the hard sphere system, complemented by the treatment of soft-core 



potentials introduced by Verlet and Weis themselves [33 ]. 

Let us emphasize that for the present purpose, approximations such as these must be re- 
garded solely as a practical and approximate mean to determine the thermodynamic property 
n£(k;n,T), which is essentially the local curvature of the free energy surface at the state 
point (n, T) 26, [m]. This property directly determines the equilibrium structure factor 
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S( e i)(k; n, T) through the equilibrium relationship nS(k; n, T)S^ eq \k; n, T) = 1, and in prac- 
tice we actually use this relationship to determine nS(k;n,T). The main message of Eq. 
(12. 3p . however, is that the experimentally observable, non-equilibrium, static structure factor 
S(k; t) is not determined by any Ornstein-Zernike equilibrium condition, but by Eq. (12. 3p 
itself, with the thermodynamic property n£(k;n,T) driving the non-equilibrium evolution 
in the manner indicated by its explicit appearance in this equation. 

B. Thermodynamic equilibrium vs. dynamically arrested states. 

In what follows, we are interested in studying the scenario revealed by the solution S(k; t) 
of Eq. ( 12. 3p . for the process of isochoric equilibration (or lack of equilibration) of the 
static structure of a system subjected to a temperature control protocol T{t) = T^W(—t) + 
T^'9(t), corresponding to a an instantaneous temperature quench to a final temperature 
denoted simply as T. Thus, the system is assumed to be prepared at an initial equilibrium 
homogeneous state characterized by a bulk particle number density n and temperature T®, 
at which its initial static structure factor is S(k; t = 0) = S^'(k). Upon suddenly changing 
the temperature of this system to the new value T, one normally expects that the system 
will reach full thermodynamic equilibrium, i.e., that the long-time asymptotic limit of S(k; t) 
will be the equilibrium static structure factor S^ eq \k]n,T) = l/n£(k;n,T). According to 
Eq. (12.31) . reaching this value is also a sufficient condition for S(k;t) to reach a stationary 
state. 

According to the same equation, however, this is not a necessary condition for the sta- 
tionarity of S(k;t), which could also be attained if lim^oo b(t) = 0, even in the absence of 
thermodynamic equilibrium (i.e., even if lim^oo S(k; t) ^ l/nS(k;n,T)). If the long-time 
stationary state attained is the thermodynamic equilibrium state, we say that the system 
is ergodic at the point (n,T). The second condition, in contrast, corresponds to dynami- 
cally arrested states, in which the long-time asymptotic limit of S(k; t) might differ from 
the expected thermodynamic equilibrium value S^ eq \k; n, T) = 1 /n£(k; n, T). Clearly, these 
are two mutually exclusive and fundamentally different classes of possible stationary states 
which can only be distinguished if we know the long-time limit of b(t). This is, however, 
not a thermodynamic property, and hence, the discrimination of the ergodic or non-ergodic 
nature of the state point (n, T) must be based on a dynamic or transport theory that allows 
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FIG. 1: Dynamic arrest line (or isodiffusivity curve with D* = 0) in the ((f), T*) state space of the 
truncated Lennard- Jones fluid, Eq. (13. ip . with v = 6. The vertical downward arrows represent two 
fundamentally different classes of irreversible isochoric processes: in the first case the fixed volume 
fraction (f> is smaller than the dynamic arrest volume fraction (j)^ (= 0.582), whereas in the second 
4> is larger than </>( a ). 

the determination of b(t). 

One such theory is precisely the SCGLE theory: to decide if the long-time stationary 
state corresponding to the point (n, T) will be an ergodic or an arrested state one can use 
the equilibrium static structure factor S^ eq \k;n,T) in Eq. ( 12. 14ft to calculate 7^ 9) (n,T). If 
the solution is infinite, we say that the asymptotic stationary state is ergodic, and hence, 
that at the point (n, T) the system will be able to reach its thermodynamic equilibrium 
state without impediment, so that lim^oo S(k; t) = \/n£{k;n,T). On the other hand, if 
the solution for ^ eq \n,T) turns out to be finite, this means that the system will become 
dynamically arrested, and that the long-time limit of S(k; t) at the point (n, T) will not 
necessarily be its thermodynamic equilibrium value S^ eq \k\n,T) = l/n£(k;n,T). Instead, 
we shall have that lim^oo S(k; t) = S^(k), with a truly non-equilibrium structure factor 
S( a '(k), different from S^ eq \k;n,T), and obtained as an alternative stationary solution of 
Eq. (12. 3p . In this manner, by calculating 7^(n, T) at all state points (n,T) one can scan 
the state space to determine the region of dynamically arrested states of the system. 

We have employed in this manner the PY/VW approximation for the equilibrium static 
structure factor S^ eq \k;n,T) of the TLJ soft-sphere model, to determine the region of its 
fluid-like ergodic states and the region of its dynamically arrested states. The resulting 
dynamic arrest transition line is represented by the solid curve in Fig. [1] for the TLJ fluid 



12 



with v — 6, whose T* — > limit coincides with the dynamic arrest volume fraction (fi^ of 
the hard sphere liquid, predicted to occur at = 0.582 |s|. As indicated in the fi gure (and 
as explained at the end of the previous section) , this transition line is also the iso-diffusivity 
curve corresponding to D* = 0. 

C. Method of solution of Ecs. (|23|) . p7H) - (p2?LOl) for equilibration. 

For concreteness, let us consider the case in which the system was initially prepared to 
be in the equilibrium state corresponding to a point (0, T*W) located in the fluid like region. 
We then have two fundamentally different possibilities, also illustrated in Fig. [U either the 
final point ((f), T*^') lies in the ergodic region of the dynamic arrest diagram, or else, it lies 
in the region of dynamically arrested states. The first case is achieved, for example, if the 
volume fraction of the isochoric irreversible process is smaller than the dynamic arrest volume 
fraction (p^ = 0.582 of the hard sphere liquid. This isochoric quench (0, T*w) — >• ((f), T*^') 
will then eventually lead to the full equilibration of the system. In the second case, in which 
the fixed volume fraction <p must be larger than = 0.582 (and the final temperature 
sufficiently low) the solution of Eqs. (12. 3ft . (I2.7p - (l2.10j) will describe the irreversible aging of 
the glass-forming liquid quenched to a point inside the dynamically arrested region. 

In either case, solving Ec. (I2.3P for S(k;t) starts with the formal solution in Eq. (12. 4p . 
written as 

S*(k; u) = S {i \k)e- a{k)u + Sf(k) (1 - e - a{k)u ) . (3.2) 

This expression interpolates S*(k;u) between its initial value S^'(k) = S^ eq '(k;(f),T*^') 
and its expected long-time equilibrium value S e j q (k) = S^ eg \k;(f),T*^) = [n^ (k)}~ x . 
Clearly, the solution S(k; t) in Eq. (12. 4 j) can be written as 

S(k;t)=S*(k;u(t)), (3.3) 

with u(t) defined in Eq. (12. 6p . The inverse function t(u) is such that u(t(u')) = u' and 
t(u(t')) = t'. The differential form of Eq. (12. 6p can be written as dt = du(t)/b(t). Upon 
integrating this equation, we have that t = J* du(t')/b(t'), which can also be written, after 
the change of the integration variable t' , to u' = u(t'), as 
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with the function b*(u) defined as b*(u) = h(t(u)). These general observations greatly 
simplify the mathematical analysis and the numerical method of solution of the full NE- 
SCGLE theory under the particular conditions considered here. 

To see this, let us consider a sequence S*(k; u n ) of snapshots of the static structure factor, 
generated by the simple expression in Eq. (13. 2p when the parameter u attains a sequence 
of equally-spaced values u n , say u n = nAu (with a prescribed Am and with n = 0,1, 2,...). 
The fact that S(k;t) can be written as S(k;t) = S*(k;u(t)) implies that this sequence will 
be identical to the sequence S(k; t n ) generated by the exact solution in Eq. (12. 4p . evaluated 
at a different sequence t n {n = 0, 1, 2, ...), i.e., at a sequence of values of the time t, given by 
t n — f n [l/b*(u')]du'. In other words, the nth member of the sequence of static structure 
factors can be labeled either with the label u n , as S*(k; u n ), or with the label t n , as S(k; t n ). 
For sufficiently small Au, the discretized form of the previous relationship between t n and 
u n can be written as 

t n +i =t n + Au/b*(u n ). (3.5) 

Thus, in practice what we do is to solve the self-consistent system of equations (I2.7[) - ( |27TTT) 
with S(k; t) replaced by each snapshot S(k; t n ) = S*(k; u n ) of the sequence of static structure 
factors. This yields, among all the other dynamic properties, the sequence of values b*(u n ) 
of the function b*(u). This sequence can then be used in the recurrence relation in Eq. 
(13. 5p to obtain the desired time sequence t n , which allows us to ascribe a well-defined time 
label to the sequence S(k; t n ) of static structure factors and to the sequence b(t n ) of the 
instantaneous mobility b(t). Of course, since the solution of equations (l2.7[) -( f27TTj) yields all 
the dynamic properties, we also have in store the corresponding sequence of snapshots of 
dynamic properties such as F(k,r;t n ), Fs(k,r;t n ), the a-relaxation time T a (t n ), etc. 

IV. EQUILIBRATION OF SOFT-SPHERE LIQUIDS. 

We have applied the protocol just described, which solves the full NE-SCGLE theory 
(Eqs. (l2.4p -f l2~TT|) ). to the description of the isochoric irreversible evolution of the structure 
and the dynamics of the TLJ soft sphere liquid, after the instantaneous quench starting 
from an equilibrium fluid state. To continue the analysis, however, it is convenient to 
discuss separately the two mutually exclusive possibilities illustrated by the two vertical 
arrows in Fig. [TJ In this section we shall concentrate on the conceptually simplest case of 
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the full equilibration of the system, and in the following section we shall discus the process 
of dynamic arrest. 

A. Ordered sequence of non-equilibrium static structure factors. 

Let us thus illustrate the isochoric quench in which both, the initial and the final points, lie 
in the fluid-like region. For concreteness, we consider a cooling process, 

y« > tV\ such that 

&W > b (f) > 0, with T* (i) = 0.1, (f) = 0.56 (< <p {a) = 0.582), and with the final temperature 
corresponding to the deepest quench, T*w) = 0. The initial and final equilibrium static 
structure factors, S®(k) = S^ eq \k; 0i, T*W) and S e /(k) = S^ eq \k; 0i,T* ( ^), are presented 
in Fig. [2j To visualize the transient non-equilibrium relaxation of S(k;t), we generate a 
sequence of snapshots S*(k;u n ) using Eq. (13. 2p with u = u n = nAu (n = 0,1,2,...) and 
with Am = O.Olfa 2 //) ] (~ l/4a(k), for k = k max , the position of the main peak of S eq {k)). 
From now on we shall use [a 2 /D°] as the time unit and a as the unit length. In Fig. [2] 
we include four representative intermediate snapshots of this sequence, corresponding to 
u/Au = 1, 3, 5, and 7. Let us emphasize that although these snapshots of the transient 
structure factor are linear combinations of two equilibrium static structure factors (namely, 
S^'(k) and Sj q (k)), they themselves represent fully non-equilibrium structures. 

Fig. [2] exhibits the fact that within the resolution Am = 0.01 employed to visualize 
S*(k; u), this non-equilibrium structure relaxes very quickly to its long-time equilibrium limit 
S eq {k) at most wave-vectors, except in two regions: in the vicinity of k max , as appreciated 
in the figure, and in the long-wavelength limit, A; — ?- 0, not apparent in the main figure, 
but illustrated and discussed below. Thus, except in these two wave-vector domains, the 
non-equilibrium snapshots of S*(k; u) shown in the figure are already indistinguishable from 
S eq {k). The fact that for large wave- vectors, k > k max , the structure S*(k; u) approaches very 
fast its final equilibrium value S eq {k) is understood by the fact that a(k) = 2k 2 D° / S e j? q (k) 
increases with k 2 while S^ik) decreases from its maximum value towards its unit value at 
large k. To the left of k max , on the other hand, although a(k) decreases with k 2 , there is a 
dramatic drop of the static structure factor from its large value at the main peak towards the 
very small value of S eq {k = 0) of a strongly incompressible liquid. In support of this proposed 
scenario, in the inset of Fig. [2]we plot the inverse relaxation constant a~ 1 (k) = Sj q (k)/2k 2 D° 
as a function of k, which clearly exhibits a dominant peak at k = k max , and a divergence 
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FIG. 2: Snapshots of the time-evolving static structure factor S*(k; u) corresponding to the quench 
0,T*«) -> (0, T*C/3) at fixed volume fraction = 0.56 and T*W = 0.1, and T*W = 0. The darker 
thick (black) solid line is the initial structure factor S*(k; u = 0) = S^\k). The lighter thick (blue) 
solid line is the asymptotic limit S*(k;u — > oo) = S eq {k). The sequence of thinner (green) solid 
lines represent S*(k;u) for u = Au, 3Au, 5Au, and 7Au, with An = 0.01 (we use [a 2 /D°] as 
the time unit). These six structure factors also correspond, according to Eq. (|3.4h . to S(k;t) for 
t = 0, 0.036, 0.1, 0.2,0.533, and oo. Inset: Inverse relaxation constant a~ l (k) = S e f q (k)/2k 2 D°. 

at k — 0. This explains the quick thermalization of S*(k;u) in both, the large wave- vector 
domain and in the moderately small wave- vector regime < k ^ k max . 

In the really small wave- vector limit k — > 0, however, the 1/k 2 divergence of u eq (k) dom- 
inates, and prevents the thermalization of S*(k;u) within finite values of u. The crossover 
from this long-wavelength perfect slowdown, to the faster moderately-small wave-vector 
regime k ^ k max , is revealed by zooming in at the small-fc behavior of the snapshots of 
S*(k;u), as illustrated in Fig. [3]^a). In contrast with this rather trivial long- wavelength 
slowing down, the slow relaxation at and around k = k max has its origin in the large value 
attained by Sj q (k ma x), i-e., in the large strength of the interparticle correlations of spatial 
extent similar to the mean distance between the particles. Thus, this slowing down of the 
main peak of the w-evolving static structure factor is a non-equilibrium manifestation of the 
so-called cage effect. In Fig. [3](b) we present a zoom of the snapshots of S*(k;u) of Fig. [2, 
exhibiting in more detail the slower relaxation of the structure at these wave-vectors. 
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FIG. 3: Zoom of the snapshots S*(k;u n ) in Fig. [2] corresponding to (a) the long- wavelength limit 
k — > and (b) to the neighborhood of the position k max of the main peak of S eq {k). 

B. Non-equilibrium u-dependence of S*(k;u) and b*(u). 

Let us notice that the simple expression for S*(k; u) in Eq. (I3.2p . which interpolates this 



This means that if we plot the static structure factor S*(k;u) as AS*(k;u) vs. the k- 
dependent variable [a(/c)ix], the results for all the wave- vectors k must collapse onto a master 
curve independent of k and of the initial and final values S^'(k) and Sj q (k). In fact, such 
a master curve will be essentially a simple exponential function. This simplicity, however, 
will be partially lost if we plot AS*(k;u) directly as a function of the parameter u, since 
such exponential function, e - "^", will decay with u at a different rate a(k) for different 
values of the wave- vector k, as illustrated in Fig. HJ^a). In fact, if we define a fc-dependent 
equilibration value u eq {k) by the condition e - a ( k ) u ^i( k ) « e ~ 5 , we have that u eg (k) = 5a _1 (/c). 
Thus, except for the arbitrary factor of 5, the inset of Fig. [2(a) exhibits the wave-vector 
dependence of u eq (k). There we see that u eq (k) attains its largest value at the wave-vector 
kmax, corresponding to the position of the main peak of S eq (k). This slowest mode im- 
poses the pace of the overall equilibration process, thus characterized by the fc-independent 
equilibration value u eq = u eq {k max ) = 5S e f 9 (k max )/2k max D°. 

The previous discussion illustrates the properties of the ordered sequence S*(k;u n ) of 
snapshots of the function S*(k;u) for equally-spaced values u n of the parameter u. This 



function of u between S^\k) and S e f q (k), may be written as 




a(k)u 



(4.1) 



17 



sequence of snapshots, however, do not fully reveal the most important features of the real 
relaxation scenario implied by the solution (12. 4p of Eq. (12. 3p . which provides S(k;t) as a 
function of the actual evolution time t. Nevertheless, since S(k;t n ) = S* (k; u{t n )) , these 
features are fully revealed by simply relabeling the referred sequence S*(k;u n ) using the 
(not equally-spaced) sequence of labels t n given by the recurrence relation in Eq. (13.51) . 
This results in the sequence S(k; t n ) of snapshots that describes the actual time evolution of 
S(k;t). In order to carry out this program, however, we must first determine the sequence 
b*(u n ) needed in the referred recurrence relation. 

As indicated before, from any sequence of snapshots S*(k;u n ), with u = nAu (n = 
0,1,2,...), we may generate a sequence b*(u n ) of values of the time-dependent mobility 
b*{u) by solving the self-consistent system of equations (I2.7p - (I2.11I) with S(k; t) replaced by 
S*(k; u n ) for each snapshot. The resulting sequence b*{u n ) is a discrete representation of the 
function b*(u), shown in Fig. IU(a), whose resolution in the parameter u may be improved 
arbitrarily by taking Au as small as needed. The first feature to notice in the result of this 
procedure is the fact that b*(u) decays monotonically from its initial value bi to its final 
value bf > 0. This implies that the system will always remain fluid-like and will have no 
impediment to reach its expected equilibrium state. We also find that the function b*(u) 
attains its asymptotic value bf for u > u eq (~ 0.2 for the quench illustrated in the figure). 

In Fig. Ufa) we also present the results for b*(u) plotted as 

Afc» S m^l. (4.2) 
Of - bi 

We see that this plot does not exhibit any simple relationship between the decay of Ab*(u) 
and the decay of AS*(k; u). In the same figure, however, the results for b*{u) are plotted as 

Plotted in this manner we observe a more apparent correlation between the decay of both, 
AS*(k max ;u) and A6*~ 1 (u), with the parameter u. This feature remains, of course, when 
these properties are expressed as functions of the actual evolution time t, as we now see. 



C. Real-time dependence of S(k;t) and b(t). 



Once we have determined the function b*(u), using the expression for t(u) in Eq. (13. 4|) . 
or its discretized version in the recursion relation of Eq. (I3.5p . we can determine the desired 
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FIG. 4: (a) The function AS*(k;u) = 1 — exp[— a(k)u] corresponding to three different wave- 
vectors, k = 6.1, 7.26 (= k max ), and 8.1, plotted as a function of the parameter u (dashed lines). 
The circles correspond to u = u eq (k) = ba~ l (k), where [1 — AS*(k; u eq {k))} = e -5 . The dotted 
and the solid dark lines represent, respectively, the functions Ab*{u) and A6*~ 1 (u) (Ecs. (|4.2p and 



^3]) ). and the solid circle corresponds to the estimated equilibration value u eq = u, 



eq \ n max 



I. The 



inset compares the function t(u) calculated according to its exact definition in Eq. (|3.4p (solid line) 
and according to the approximation in Eq. ()4.6p . (b) Same information as in (a), but now plotted 
as a function of the actual evolution time t = b ,^ u ,^ du' . In the main figure we plot AS(k; t) and 
in the inset we plot S(k;t). 

real-time evolution of S(k;t) and b(t). In this manner we determine that in our illustrative 
example the sequence u = 0,0.01, 0.03, 0.05, and 0.07 corresponds to the sequence t = 
0, 0.036, 0.1, 0.2, and 0.533. In the inset of Fig. 11(b) we present the resulting time-evolution 
of S(k; t) for the same three wave- vectors as in Fig. E(a). This inset emphasizes the fact that 
S(k; t) evolves monotonically from its initial value S^(k) to its final value Sj q (k), sometimes 
increasing and sometimes decreasing, depending on the wave-vector considered. In order to 
exhibit a less detail- dependent scenario, in the main frame of Fig. 11(b) we present the same 
information, but formatted as AS(k;t), which is the re-labeled version (u — > t — t(u)) of 
AS*(k;u) in Eq. f !4.ip . namely, as 

S(k;t) - S®(k) 



AS(k;t) = 



(4.4) 



We similarly relabel the definitions of Ab*(u) and A6* _1 (m) in Eqs. (14. 2 p and (14. 3 p to define 
the functions Ab(t) and A6 _1 (t), which are also plotted in Fig. 11(b). The comparison of this 
figure with Fig. H(a) indicates that, except for the stretched metric of t, the overall scenario 
described by the w-dependence illustrated in Fig. @|a) is preserved in the t-dependence 
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illustrated in Fig. @|b). 

Let us notice in particular that the existence of the equilibration value u eq of the parameter 
u, beyond which b*(u) ~ 6/, allows us to define an equilibration time, t eq , as the time that 
corresponds to u eq through Eq. ( 13. 4p . 

For our specific illustrative example, this yields t eq w 48. The fact that b*(u) ~ bf for 
u ^ u eq implies, according to Eq. (13.4p . that for u > u eq the function t(u) will be linear in 
u, i.e., 

*(u) « -a{u eq ) + fej 1 !*, (4.6) 

with a{u eq ) = Jq [1/6/ — 1/6* (it')]cf?/. In the inset of Fig. @Ja) we compare this asymptotic 
expression, applied to our illustrative case (for which bj 1 = 285 and a(u eq ) = 12.8), with 
the actual t(u) calculated from Eq. ( 13.41) . 



D. Irreversibly-evolving dynamics. 

Since for each snapshot of the static structure factor S(k; t) the solution of Eqs. (12. 7ft - 
(12. lip determines a snapshot of each of the dynamic properties of the system, the process of 
equilibration may also be observed, for example, in the t-evolution of the collective and self 
intermediate scattering functions, F(k,r;t) and F s (k,r;t). In Fig. [S^a) we illustrate this 
irreversible time-evolution with the snapshots of the self-ISF Fs(k max ,T;t), corresponding 
to the same set of evolution times t n as the snapshots of S(k] t) in Fig. El We see that 
the function Fs(k,r;t) starts from its initial value Fg(k,r;t = 0) = Fg q {k, r; 0, Tj = 0.1), 
and quickly evolves with waiting time t towards the vicinity of its final equilibrium value 
F s (k } T;t = oo) = F e s q (k } r; 0,T/ = O). 

The equilibration process of Fs(k,r a ]t) can be best summarized in terms of the depen- 
dence of the a- relaxation time T a (k; t) as a function of the evolution time t. The a-relaxation 
time may be defined by the condition 

F s (k,T a ,t) = l/e. (4.7) 

The dependence of T a (k;t) on the evolution time t can be extracted from a sequence of 
snapshots of F s (k, r a ; t), such as those in Fig. E{a). The results are illustrated in Fig. [5(b), 
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FIG. 5: (a) Snapshots of the non-equilibrium self intermediate scattering function Fs(k,T;t) 
at k = k max , corresponding to the equilibration process in Fig. [2] for evolution times t = 
0, 0.036, 0.1, 0.2,0.533 and oo. (b) Non-equilibrium evolution of the dimensionless a-relaxation 
time, displayed as r a {k;t) itself (inset) and formatted as in Eq. 14.81 (main figure). The solid circle 
represents the equilibration point (t eq ,T a (k max ;t eq )). 

in which the solid line corresponds to T a (k max ;t). The solid circle indicates the crossover 
from the t-regime where T a (k max ;t) is still in the process of equilibration, to the regime 
where it has reached its final equilibrium value Tq (k max ) = T^ q (k max ] <f>,Tf). In the inset of 
the figure we plot T a (k max ; t) itself and in the main figure we plot the same information, but 
formatted as 

r a (k;t)-r^(k) 



Ar Q (fc;t) 



(4.8) 



T<P(k)-Tg\ky 

with T$(k) = r^i^&Ti). In the same figure we also exhibit similar results corresponding 
to two additional wave- vectors, different from k max (dashed lines). These results show that 
the equilibration time of T a (k;t) for these three wave-vectors is largely independent of k, 
and can be well approximated by the equilibration time t eq defined in Eq. (14. 5p , in contrast 
with the notorious wave- vector dependence of the predicted evolution of the static structure 
factor illustrated in Fig. 13(b). 

Let us finally mention another theoretical prediction regarding the kinetics of the equi- 
libration process. This refers to the similarity of the equilibration kinetics exhibited by 
the time-dependent mobility b(t), the a-relaxation time r a (k; t) at all wave- vectors, and the 
static structure factor S(k max ;t) at the wave- vector k max , when plotted in terms of the re- 
duced properties Ab(t), Ar a (k; t), and AS(k max ; t). This similarity is exhibited in Fig. |6]for 
the illustrative quench at fixed <fi = 0.56, and means that indeed the evolution of S(k max ; t) 
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FIG. 6: Time-evolving a-relaxation time T a (k max ;t) (solid line) mobility b(t) (dash line) and 
static structure factor S(k max ;t) (dotted line) as a function of evolution time t, and expressed as 
AT a (k max ;t), A6 _1 (t), and AS(k max ;t). 

(which is slower than the evolution of S(k; t) for other wave-vectors) sets the overall relax- 
ation rate exhibited by the dynamic properties Ab(t) and Ar a (k; t). Thus, from this point of 
view, we may use either of these characteristic dynamic properties to describe the predicted 
kinetics of the equilibration process. 

E. Dependence on the initial temperature of the quench. 

Up to this point we have illustrated the main features of the isochoric quench (0, T*W) — > 
(0, r*w) at fixed volume fraction = 0.56, using for concreteness the values T*W = o.l 
and T*tf) = 0. We are now ready to analyze how the scenario just described depends 
on the initial temperature T*W and on the volume fraction at which we perform the 
quench. Let us start by considering the dependence on T*M. Rather than attempting a 
comprehensive illustration of this dependence in terms of the evolution of the static structure 
factor S(k; t) and of the various dynamic properties, we use the dimensionless mobility b(t) as 
a representative property bearing the essential information about the equilibration process. 
This fc-independent property determines the mapping from the parameter u to the real time 
t, through the definition of the functions u(t) and t{u) in Eqs. (12. 61) and (13 .41) . 

Thus, in Fig. [7(a) we present plots of &~ x (t) as a function of t for three representative 
values of the initial temperature T*W, namely, T*W= 1.0, 0.1, and 0.01, keeping the same 
final temperature T*^' = and the same volume fraction 0=0.56. This figure reveals 
two remarkable features. In the first place, the equilibration time t eq seems to be rather 
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FIG. 7: Non-equilibrium, time-dependent (a) mobility b(t) and (b) a-relaxation time T a (t), as a 
function of evolution time t, for the isochoric quench (</>, T*w) — T*(^ = 0) at fixed volume 
fraction ^ = 0.56, for three different initial temperatures, T*W= 1.0, 0.1, and 0.01. The insets 
exhibit a power law fit of the transient before saturation. 

insensitive to the temperature J^W of the initial state. In other words, the system will 
reach the final equilibrium state in about the same time, t eq w 48, no matter if the initial 
temperature is T*W= l.o, 0.1, or 0.01. To emphasize this feature we have highlighted the 
common equilibration point of the three curves. The second remarkable feature is that 
during the transient stage of the equilibration process, the evolution of as a function 

of t follows approximately a power law, ~ At x , with the exponent x and the amplitude 

A depending on the initial temperature T*M. i n the inset of the figure we exhibit the power 
law fit of the transient, indicating the resulting value of the exponent x and amplitude A. 

Exactly the same trend is also reflected in the evolution of the intermediate scattering 
function, as observed in the results for the a-relaxation time shown in Fig. [3(b). This 
information is important, since many times it is this dynamic parameter what is monitored 
in simulations and in some experiments. 

F. Dependence on the volume fraction of the quench. 

Let us now discuss the dependence of the equilibration process on the volume fraction <p. 
Once again we first use the time evolution of b(t) to illustrate this dependence. In Fig. E^a) 
we plot b~ l (i) as a function of t for a set of values of the volume fraction 0, corresponding 
to the metastable regime of the hard-sphere liquid. According to these results, the inverse 
mobility 6 _1 (£;</>) reaches its equilibrium value b"7 (<f>) after a (^-dependent equilibration 
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time t eq (<f)). To emphasize this prediction, the solid circles in the figure highlight the points 
(t eq ((f)) , bj 1 ((p)) . These highlighted points, as indicated in the figure, align themselves to a 
good approximation along the dashed line of the figure, corresponding to the approximate 
relationship bj 1 ^) « 4 x [t e « (</>)] L05 

This relationship between the equilibration time t eq ((f>) and bf((f>) is one of the most 
remarkable predictions of the present theory, bearing profound physical implications. To 
see this let us recall that the dimensionless mobility bf((f>) is just the scaled long-time self- 
diffusion coefficient D*((p,Tf) = D L ((f>,Tf)/ D° of the fully equilibrated system at the final 
point ((f), Tf), which for the present isochoric quench down to zero temperature, Tf = 0, is 
the dimensionless equilibrium long-time self- diffusion coefficient of the hard sphere liquid, 
D* HS ((f>) = D*((f>,Tf = 0). This property can be calculated using the equilibrium version 
of the present theory [22j and, as discussed below (see fig. E^b)), such calculation leads to 
the prediction that D* HS (<f)) vanishes at (f>^ = 0.582, according to the power law D* HS ((f)) oc 
((p^ - (f>) 2 - 2 . As a consequence, if t eq ((f>) « 0.25 x bj 1 ((f>) (oc D* H g((f>)), we must expect that 
as (p (p^ the equilibration time will diverge according to t eq ((f)) oc ((f)^ — (f))~ 2 ' 2 - 

This predicted divergence of the equilibration time constitutes a strong and interesting 
proposal, which requires, of course, a critical assessment and validation. We shall return to 
this discussion later on in the paper, but at this point, let us carry out a similar analysis, now 
using the a-relaxation time r a (t;0) (whenever we omit the wave- vector k as argument of 
r a (k, t; (f)) is because a specific value for k is being assumed fixed, most frequently k ~ k max ). 
Thus, in Fig. E](b) we plot r a (t] (f>) as a function of t for the same set of values of the volume 
fraction (f> as in Fig. El^a). Here again we find that T a (t](f>) reaches its equilibrium value 
T^ q ((f>) after the same evolution time t eq ((p) as in the case of b(t;(f>). Also here the solid 
circles highlight the points (t eq ((f>) , r^ q ((f>)) , and the dashed line of the figure imply that 
t eq ((f>) ~ 40 x [Ta 9 ] o ' 95 (0)- This implies that the time t eq ((f>) required to equilibrate the 
system will grow at least about as fast as the equilibrium value r^ q of the a-relaxation time, 
and that both properties increase strongly with (f>. 

Although one can discuss additional features of the class of irreversible process corre- 
sponding to the full isochoric equilibration of the system after its sudden cooling, it is now 
important to contrast the scenario just described, with that of the second class of irreversible 
processes. This involves the dynamic arrest of the system, and is the subject of the following 
section. 
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FIG. 8: Non-equilibrium, time-dependent (a) mobility and (b) a-relaxation time T a (t), as a 

function of evolution time t for various volume fractions corresponding to the metastable regime 
of the hard-sphere liquid. The arrows indicate the direction of increasing <f>. The dashed lines are 
the power laws bj 1 ^) «4x [i e 9(0)] L05 and T% q {<t>) « 0.025 x [t e <?(»] L05 . 

V. AGING OF SOFT-SPHERE LIQUIDS. 

Let us recall at this point that the NE-SCGLE description of the spontaneous evolution 
of the structure and dynamics of an instantaneously and homogeneously quenched liquid is 
provided by the simultaneous solution of Ecs. f l2.3p - fl2.10p . As discussed in subsection IIIIBl 
there exist two fundamentally different classes of irreversible isochoric processes, represented 
by the vertical downward arrows in Fig. [TJ In the previous section we described the resulting 
scenario for the most familiar of them, namely, the full isochoric equilibration of the system. 
In this section we present the NE-SCGLE description of the second class of irreversible 
isochoric processes, in which the system starts in an ergodic state and ends in the region 
where it is expected to become dynamically arrested. 

Thus, let us continue considering the TLJ model system introduced in Sect. HH] (Eq. 13.11 
with v = 6), subjected to the sudden isochoric cooling at fixed volume fraction <fi = 0.6, 
larger than <j>^ = 0.582, from the point (</>, T*W = 0.1) in the ergodic region, to the point 
(0,T*W = 0) in the region of dynamically arrested states. By construction, the solution 
7W of Eq. fl2.14p . obtained using S^ eg \k; 0, T*^>) as the structural input, is 7W = 00. In 
this sense, the present class of process is identical to the first one, discussed in the previous 
section. The main difference lies, of course, in the fact that in the present case the solution 
of Eq. (I2.14p for the squared localization length 7^, obtained using the equilibrium static 
structure factor S^ eq \k; (f),T*u>) of the final point as input, will now have a finite value. 
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To see the consequences of this difference, let us go back to subsection IIII C] and consider 
the function S*(k; u) in Eq. (13. 2p . with < u < oo. For each value of u we may use S*(k; u) 
in the bifurcation equation (I2.14p for 7(i), now denoted as 7* (it)- Throughout the previous 
section it was implicitly assumed that 7* (it) = 00 for < u < 00, an assumption based 
on the fact that the system started and ended in a fluid-like state. In the present case, 
however, although the system starts with the condition that 7*(w = 0) = 00, we know that 
the final point ((f), T*^ = 0) corresponds to an arrested state, so that 7^ = 7*(u = 00) has 
a finite value. This means that somewhere between u = and u = 00 the function 7*(w) 
changed from infinity to a finite value, and this then implies the existence of a finite value 
of u, such that 7*(w) remains infinite only within the interval < u < u^ a \ Thus, in 
the present case the simultaneous solution of Ecs. (I2.3l) - (l2.10p starts in practice with the 
precise determination of vS a \ 

A. Method of solution of Ecs. (HO) . (l277ft - d2?L0ft for aging. 

To determine the critical value vS a \ let us consider again the sequence S*(k;u n ) of snap- 
shots of the static structure factor, generated by the expression in Eq. (I3.2p with u n = nAu 
(n = 0, 1,2, ...). Since we have assumed that initially the system is fluid-like, the value of 
cannot be u^ a ' = 0. Thus, let us employ each snapshot of the sequence S*(k;u n ), with 
n — 1,2, as the static input of Eq. (12.141) . thus determining the sequence 7* = 7*(u n ) 
of values of 7*(w), which starts with 7q = 00. If 7* turns out to be finite, then one may 
take a smaller w-step Am, until this does not happen. For a sufficiently small Au, there 
will be an integer n a such that 7* = 00 for n < n a and 7* is finite for n > n a , i.e., such 
that u na < < W(„ a+ x). This process can be refined by decreasing Au, so that one can 
determine with arbitrary precision for the given initial and final conditions ((f), T*W) and 
(0,T*( / )). For example, one can readily perform this procedure for the quench indicated by 
the right arrow of Fig. [1] (from the point ((f), T*^ = 0.1) to the final point ((f), T*^ = 0) at 
fixed = 0.6), with the result w (a) = 0.0128. 

Once one has determined with the desired precision, one can construct a new sequence 
ui of (N + 1) equally-spaced values of u, defined as ui = / x (vS a ' /N) with < I < N, along 
with the corresponding sequence S*(k;ui) of snapshots of the static structure factor (using 
Eq. (13. 2p ). Since the sequence S*(k;ui) is identical to the sequence S(k;ti) (with t\ such 
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FIG. 9: (a) Mobility function b*(u) plotted as a function of the parameter u and of the difference 
_ u ) (inset) for the sudden isochoric cooling from the point (4> = 0.6, T*W = 0.1) to the point 
O = 0.6, T*(/) = 0). (b) Scaled long-time self-diffusion coefficient D* HS (<p) = D L {4))/D of the 
hard-sphere liquid as a function of volume fraction <p and of the difference (cj)^ — 4>) (inset). 



that U[ = b(t')dt'), to each member of this sequence, the self-consistent Eqs. fl2.8|) - fl2.10p 
assigns a snapshot of the full dynamics of the system. In particular, the use of Eq. (12. 7p 
generates a sequence of values b*(ui) of the mobility b*(u), with arbitrary resolution (set by 
the number iV of w-steps). To illustrate these concepts, in Fig. [9|a) we present the results 
for b*{u) corresponding to the specific quench under discussion. Notice that, as expected, 
b*(u) — > as u approaches from below. 

A simple ansatz to model this limiting behavior is 

b*(u) ^B (u {a) -uf. (5.1) 

In the inset of Fig. [9](a) we plot b*(u) vs. (u^ — u) to determine the value of the exponent 
j3 and the pre-factor B Q , with the result /3= 2.2 and B = 9.5. We performed similar calcu- 
lations varying the initial temperature T*W, and found the value (3= 2.2 of the exponent is 
independent of T*W, so that the dependence of b*(u) on the initial temperature is carried 
only in the pre-factor B . For example, we found that B (T*^) = 9.5, 33, and 490, for 
T*«= 0.1, 0.05, and 0.01, respectively. 

Another remarkable feature of the w-dependence of b*(u) illustrated in the inset of Fig. 
Efa) is its similarity with the volume fraction dependence of the scaled long-time self-diffusion 
coefficient of the fully equilibrated hard-sphere system, D* HS ((j)) = D L (<fi,Tf = 0)/D°. This 



property can be calculated using the equilibrium version of the SCGLE theory [22j, and 
the results are exhibited in Fig. 0(b). As discussed before 
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231 ]. the theoretical prediction 
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FIG. 10: Non-equilibrium mobility b(t) as a function of evolution time t for the quench processes 
from the point ((f) = 0.6, T*® = 0.1) to the point (<j) = 0.6, T*W = 0). The inset exhibits the long-t 
asymptotic decay b(t) ~ bot~ v described by Eq. (|5.6|) . 

is that D* HS ((p) vanishes at the dynamic arrest volume fraction = 0.582. The results 
of Fig. EEb) show that in the vicinity of <jt- a \ the function D* HS (cf)) follows the power law 
D* HS (cf)) oc (4>^ — 0) 2 ' 2 , i.e., it vanishes at with the same exponent as b*{u) vanishes at 
u = u( a \ 

1. Asymptotic decay b(t) oc t~ v . 

At this point let us notice that the sequence b*(ui) must be identical to the sequence 
b(ti) of values of b(t) at the times U = f^ l [l/b*(u')]du'. The sequence of times U can be 
determined by means of the approximate recurrence relationship in Eq. (13.51) . i.e., 

t i+1 = t l + Au/b*(u l ), 

with Au = u^'/N, and it allows us to transform the sequence b*(u{) into the discrete 
representation b(t{) of the function b(t). The results for b(t) are plotted in Fig. [TQl to 
exhibit the fundamentally different behavior of the functions b*(u) and b(t). While the 
former has a well-defined zero at a finite value of its argument, namely, at u — u^ a \ the 
function b(t) decays to zero in a much slower fashion. In fact, as we now discuss, one of the 
main predictions of the NE-SCGLE theory is that b(t) will remain finite for any finite time 
t, and only at t = oo the mobility will reach its asymptotic value of zero. Thus, the system 
in principle will always remain fluid-like, and the dynamic arrest condition b(t) = will only 
be reached after an infinite waiting time. 
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Let us actually demonstrate that the value of t corresponding to is = oo, and 
that the mobility b(t) decays as a power law with t. To discuss the first issue, let us recall 
Eq. ( 13. 4p . which writes the function u(t) as 

t(u) = L ¥WY (5 ' 2) 

where the function b*(u) is, of course, b*(u) = b(t(u)). According to this result, and to Eq. 
( 15. ip . we can write 

for u in some vicinity uq < u < of u^. This implies that, if the exponent /3 is larger 
than unity, then t(u) will diverge as u approaches according to 

^ (g-ik ' (5 ' 4) 

As a consequence, the dynamic arrest time = t(u^ a ') will be infinite, which is what we 
set out to demonstrate. 

Let us now discuss the possibility that b(t) decays as a power law with t. For this, let us 
invert the function t(u) in the previous equation, and write it as 

u(t) ^u {a) - {(/3 -l)B ty^T). (5.5) 

Since, according to Eq. (12. 6p . b(t) = du(t)/dt, the time derivative of this asymptotic expres- 
sion will yield the asymptotic form for b(t), namely, 

b(t) « b r v , (5.6) 

with 

6 = p-l)^o]-^, (5.7) 

and 

" = jfh) ( OT <"-!) = (^T))' < 5 - 8 > 

The latter result implies that if one of the exponents (/3 or if) is larger than unity, then 
the other is also larger than unity. It also implies that if one of them is larger than 2, then 
the other is smaller than 2, and viceversa. In the inset of Fig. [10] we compare the actual 
NE-SCGLE results for b(t) in the main figure, with the approximate asymptotic expression 
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in Eq. (15. 6 j) with a fitted exponent 77, with the result that 77 = 1.83. This value coincides 
with the expected result 77 = /3/(/3 — 1) with /3 = 2.2. As indicated above, we performed 
similar calculations varying the initial temperature T*W ; and found that the scenario just 
described is indeed independent of T*W. Thus, in the asymptotic expression in Eq. (15. 6p 
only the pre-f actor 6 depends on T*W ; anc l the approximate expression in Eq. (15.71) provides 
an indicative estimate of its actual value. 

2. Dynamically arrested evolution of S(k;t). 

The properties of the non-equilibrium mobility function b(t) that we have just described 
reveals the main feature of the time evolution of the static structure factor S(k; t) when the 
system is driven to a point (<f),T*) in the region of dynamically arrested states. We refer to 
the fact that under such conditions, the long-time asymptotic limit of S(k; t) will no longer 
be the expected equilibrium static structure factor S^ eq \k; 0, T*), but another, well-defined 
non-equilibrium static structure factor S^ a \k), given by 

S( a \k) = S (i \k)e- a(k)u(a) + S e /(k) (l - e-^ 1 "'"') . (5.9) 

This non-equilibrium static structure factor not only depends on the final point (0, T*), but 
also on the protocol of the quench (in the present instantaneous isochoric quench, this means 
on the initial temperature T*w). 

To see the emergence of this scenario, let us consider the sequence S*(k; u{) of snapshots 
of the static structure factor generated with Eq. (13. 2p . for the finite sequence u\ of (N + 1) 
equally-spaced values of u defined as ui = I x (u^ a '/N) with < I < N. According to 
Eq. (15.21) . and to its asymptotic version in Eq. (15.41) . in the present case the finite range 
< u < maps onto the infinite physically relevant range < t < 00 of the evolution 
time t (in contrast with the equilibration processes studied in the previous section, in which 
the infinite range < u < 00 maps onto the infinite range < t < 00). Since the 
sequence S*(k;ui) is identical to the sequence S(k;ti), with t\ = J^ 1 du' /b*(u'), then the 
sequence of snapshots S(k; ti) describing the full evolution of S(k; t) will be generated by 
a sequence of snapshots of S*(k;u) with u only in the range < u < u^ a \ In other 
words, in the present case none of the snapshots of S*(k;u) with u > will map onto 
any physically observable snapshot of S(k] t), and this applies in particular to the snapshot 
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S*(k] u = oo), corresponding to the expected equilibrium static structure factor S^ eq \k). In 
this manner, the long time limit of S(k;t), normally being the ordinary equilibrium value 
S^ eq \k) (= \nE^\k)]~ l ) , is now replaced by a non-equilibrium dynamically arrested static 
structure factor S^(k) given, according to Eq. (I3.2p . by the expression in Eq. (I5.9p . 

Besides the remarkable prediction of the existence of this well-defined non-equilibrium 
asymptotic limit of S(k;t), the second relevant feature refers to the kinetics of S(k;t) as it 
approaches S^ a \k). To exhibit this feature, let us subtract Eq. ( 15. 9 p from Eq. (I2.4j) . This 
leads to 

S(k; t) - S {a \k) = A(k) \ e -^)Ht)-u^] _ J ^ ( 5 _ 1Q ) 

with 

A(k) = e- a{k)u(a) {S {l \k) - S e /(k)}. (5.11) 
At long times, when [u(t) — u^ a '] is small, this equation reads 

S(k;t) -S {a \k) « A{k)a{k)[u (a) - u{t)}. (5.12) 

From Eq. (15. 5p . however, we have that u {a) -u(t) sa {{P-l)B Q t} c 3 - 1 ), so that the previous 
long-time expression for S(k; t) can be written as 

S(k; t) - S (a \k) w D(k)t~Ti£v, (5.13) 

with 

D(k) = A(k)a(k){(P - l)B }~v£v . (5.14) 

Thus, we conclude that, contrary to the kinetics of the equilibration process, in which 
S(k;t) approaches S^ eq \k) in an exponential-like fashion, this time the decay of S(k;t) to 
its stationary value S^ a '(k) follows a power law. At very short times, however, b(t) ~ 
and hence, u(t) ~ b^'t. Thus, according to Eqs. (12 .4p and (12.61) . we have that the very 
initial evolution of S(k; t) might seem to approach its expected equilibrium value S^ eq '(k) = 
\nE^\k)\~ l in an apparently "exponential" manner, with a relaxation time t app ~ l/a(k)b^\ 
This apparent initial exponential evolution, however, crosses over very soon to the much 
slower long-time evolution of S(k;t) described by the asymptotic expression in Eq. (15.141) . 

Fig. CD] illustrates with a sequence of snapshots the predicted non-equilibrium evolution 
of S(k;t) after the isochoric quench at = 0.6 from T*W = 0.1 to T*v> = 0. There we 
highlight the initial static structure factor S^(k) = S^ eq \k;(p,T*^) and the dynamically 
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FIG. 11: Snapshots of the non-equilibrium evolution of S(k;t) (thin solid (green) lines) correspond- 
ing to the isochoric quench at fixed volume fraction <f> = 0.6, from T*W = 0.1 to T*w) = 0. The 
dashed (black) line is the initial structure factor S^(k). The dotted (blue) line is S e ^(k). The 
thick solid (red) line is the dynamically arrested asymptotic solution S a (k) of Eq. (|2.3p . given by 
Eq. (|5.9p . In the inset, the solid line is the maximum of S(k; t) as a function of the evolution time 
t, and the dashed line is the maximum of [S« (k)e~ a W b{i)t + [nS^^k)}- 1 ^ - e -°( fc ) 6(l) *)]. 

arrested long-time asymptotic limit S^ a \k) of the non-equilibrium evolution of S(k;t). For 
reference, we also plot the expected, but inaccessible, equilibrium static structure factor 
S^\k] </>,T*W)) = l/n£(k;(j),T*W) corresponding to the final temperature T*W = 0. Re- 
garding the kinetics of the non-equilibrium evolution, in the inset we plot the evolution of 
the maximum of S(k; t) as a function of t to illustrate the fact that S(k; t) approaches S^ a \k) 
much more slowly, in fact as the power law [S^(k) — S(k;t)] oc t -a83 . For reference, we 
also plot the maximum of the function [S^(k)e- a( - k ^ h + [^^(A;)]- 1 ^ - e - Q W feW *)] which, 
according to Eq. ( 13. 21) . would describe the evolution of S(k;t) if b(t) remained constant, 
b(t) = 6W. 

B. Aging of the dynamics 

Let us now discuss how the scenario just described manifests itself in the non-equilibrium 
evolution of the dynamics. We first recall that for each snapshot of the static structure 
factor S(k;t), the solution of Eqs. fl2.7p - fl2.lip determines a snapshot at waiting time t of 
each of the dynamic properties of the system. Thus, the process of dynamic arrest may 
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also be observed, for example, in terms of the t-evolution of the self intermediate scattering 
function F$(k, r; t) or of the a- relaxation time T a (k; t). In Fig. [T27 a) we present a sequence of 
snapshots of the ISF F$(k, r; t) (thin solid lines), evaluated at the fixed wave- vector k = 7.1, 
plotted as a function of correlation time r, for a sequence of waiting times t after the sudden 
temperature quench from T*W = 0.1 to T*^ = at fixed volume fraction = 0.6. 

In the figure we highlight with the dashed line the initial ISF Fs(k,r;t = 0). The 
(arrested) non-equilibrium asymptotic limit Fg a \k,r) = lim^oo F s (k, r; t) is indicated by 
the solid line, whereas the dotted line denotes the inaccessible equilibrium ISF Fg 29 \k } T) } 
i.e., the solution of Eqs. f l2.8l) -( l2TTTi) in which the final equilibrium static structure factor 
S^ eq \k]4> = 0.6, T*W = 0) (also inaccessible) is employed as static input. We observe that 
at t — 0, Fs(k,T;t) shows no trace of dynamic arrest, but as the waiting time t increases, 
its relaxation time increases as well. In the figure we had to stop at a finite waiting time, 
but the theory predicts that the ISF Fs(k,r;t) will always decay to zero for any finite 
waiting time t, and continues to evolve forever, yielding always a finite, ever-increasing, a- 
relaxation time r a (k; t). The relaxation of F$(k, r; t) is characterized by a fast initial decay 
(/?- relaxation) to an increasingly better defined plateau, whose height fo{k) is not determined 
by the expected equilibrium ISF Fg q (k,r), but by the non-equilibrium asymptotic limit 
Fg a \k,r). In other words, fo(k) is the "true" non-equilibrium non-ergodicity parameter 
/^(^^lim^F^^r). 

From this sequence of snapshots of Fs(k,r a ;t) we can extract the t-evolution of the a- 
relaxation time T a (k; t) defined in Eq. (14. 7p . The results allows us to notice one of the main 
features of the predicted long-r decay of Fs(k,r;t), namely, the long-time collapse of the 
curves representing F$(k, r; t), corresponding to different evolution times t (like those in Fig. 
rnZT a)). onto the same stretched-exponential curve upon scaling the correlation time r with 
the corresponding T a (k; t). In other words, at long times F s (k, r; t) scales as 

F s (k,r;t) w he~<^))\ (5.15) 

where fo is the height of the plateau of Fg a \k, r) and a = 1 + In f (so that Fs(k } r a ; t) = 
e _1 ), and with (3 being a fitting parameter. This scaling is illustrated in Fig. fT2T b) with 
the sequence of results for Fs(k, r; t) in Fig. fT2T a) now plotted in this scaled manner, which 
are then well represented by the stretched-exponential function above, with fo = 0.624, 
a Q = 0.528 and = 0.82. 
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FIG. 12: (a) Sequence of snapshots of the intermediate scattering function Fs(k,T;t) at k = 7.1 
(thin solid (green) lines) plotted as a function of correlation time r for a sequence of values of 
the waiting time t (= 0.25, 5.6, 106, 400, 590, and 1600) after the sudden temperature quench 
at fixed volume fraction (ft = 0.6, from T*w = 0.1 to T*"> = 0. The dashed line is the initial 
ISF Fs(k,r;t = 0), the thick solid (red) line is the non-equilibrium asymptotic limit Fg\k,r) = 
Fs(k,T;t = oo), and the dotted (blue) line is the expected (but inaccessible) equilibrium ISF 
corresponding to the hard-sphere system at (ft = 0.6. (b) Same sequence of snapshots of Fs(k, r; t) 
plotted as a function of the time t scaled with the a-relaxation time r a (t). The dot-dashed line 
here is the stretched exponential 0.624 x exp [— 0.528(t/r Q ) a82 ] 

The non-equilibrium evolution of the dynamics can be summarized by plotting T a (k; t) as a 
function of waiting time t. This is done here in Fig. Q31 where we plot r a (t) (= r a (k = 7.1, t)) 
as a function of t. The thick dark solid line in Fig. [T3Ta) and (b) derive from the sequence 
of snapshots of F s (k, r a ; t) in Fig. [121(a). corresponding to the quench at (ft = 0.6 with initial 
temperature T"*W = 0.1. As indicated in these figures, at long waiting times we find that 
r a (k; t) increases with t according to a power law that is numerically indistinguishable from 
r a (t) oc f with i] m 1.83. In other words the present theory predicts, taking into account 
Eq. (15.61) . that at long waiting times, r a (k;t) diverges with t with the same power law as 
b-\t). 

Besides these results, Fig. [TST a) also presents theoretical results for two additional quench 
programs that differ only in the initial temperature, namely, T*w = 0.05 and J^w = 0.01. 
These three initial temperatures lie above the dynamic arrest transition temperature T*^ ((ft) 
corresponding to the isochore (ft = 0.6, which is T*^ a \(ft = 0.6) = 0.004. The first feature 
to notice is that the detailed waiting time dependence of r a at short times may be strongly 
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FIG. 13: Waiting-time dependence of the a-relaxation time T a (t; J'*W ; 0) (defined in the text) after 
the sudden temperature quench at fixed volume fraction cj>, from an initial temperature T*W to a 
final temperature T*w) = 0. In (a) we present the results for the initial temperatures T*W =0.1, 
0.05, and 0.001 at the same volume fraction <fi = 0.6. In (b) we fix the initial temperature T*W =0.1 
and present results for <j) = 0.6 and other volume fractions. The dashed lines indicate the asymptotic 
power-law r a {t) oc At x that fits the results in the indicated regimes. 



quench- dependent, but the asymptotic power law r a (t) oc t v with the exponent r) « 1.83 is 
independent of the initial temperature T*M. Complementing this information, Fig. [T3T b) 
describes the dependence of the evolution of r a (t; T*^\ <p) on the value of the volume fraction 
at which these isochoric processes occur, assuming that each of them start and end at the 
same initial and final temperatures, T*M = 0.1 and T*w = 0. The main feature to notice in 
these results is that the long-time asymptotic growth of r a (t; <fi) with waiting time t is also 
characterized by the power law T a (t) oc t L83 . 



VI. CROSSOVER FROM EQUILIBRATION TO AGING. 

Of course, one could continue describing the predictions of the NE-SCGLE theory re- 
garding the detailed evolution of each relevant structural and dynamic property of the 
glass-forming system along the process of equilibration or aging. At this point, however, 
we would like to unite the main results of the previous two sections in a single integrated 
scenario that provides a more vivid physical picture of the predictions of the present theory. 
With this intention, in Fig. HUf a) we have put together the results for the isochoric evolution 
of T a (t; (ft) previously presented in Figs. Mjo) and HBT b). corresponding to the quench at fixed 
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FIG. 14: (a) Waiting time dependence of the a-relaxation time T a (t;4>) for a sequence of fixed 
volume fractions (the arrow indicates increasing <p). (b) ^-dependence of the a-relaxation time 
T a (t;4>) for the sequence of fixed waiting times t = 10°, 10 1 , 10 2 , 10 3 , and 10 4 (from bottom to 
top). The dashed line is t^{4>) = lim^oo T a (t; (p), which is the equilibrium a-relaxation time of 
the hard-sphere system, predicted by the equilibrium SCGLE theory. The solid circles highlight 
the crossover points (4>^ c \t), Ta q (<fi^)) for each waiting time t shown. The inset of (b) shows the 
evolution of the crossover volume fraction dy& (t) predicted by the NE-SCGLE theory (solid line) 
and determined in the simulations of Ref. [8| (solid squares) . 

volume fraction 0, from an initial temperature T*^> = 0.1 to a final temperature T*v> = 0, 
for volume fractions (f> smaller and larger than (jy- a > = 0.582. 

Displaying together these results allows us to have a richer and more comprehensive sce- 
nario of the transition from equilibration processes to aging processes in the soft-sphere 
glass-forming liquid, discussed separately in the previous two sections. According to the 
NE-SCGLE theory, the dynamic arrest transition is in principle a discontinuous transition, 
involving the abrupt passage from one pattern of evolution (equilibration) to the other (ag- 
ing) when the control parameter <p crosses the singular value (p^ = 0.582. The discontinuous 
nature of this kinetic transition is rooted in the abrupt transition contained in the equilib- 
rium version of the SCGLE theory, which actually predicts the existence and location of the 
dynamic arrest transition line (see Fig. [T3|) . The zero-temperature limit of this transition 
line corresponds to the critical volume fraction (j>^ = 0.582. Thus, the evolution of the 
system after the temperature quench from an initial temperature J"*( l ) = 0.1 to a final tem- 
perature T*U) = is dramatically different if the volume fraction of the isochoric process is 
smaller or larger than this critical volume fraction. 

However, in order to actually witness this dramatic difference, we would have to perform 
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observations at volume fractions infinitesimally closer to <fi( a \ and within an evolution time 
window much larger (in fact, infinite) than in any real experiment or simulation. In fact, 
what we would like to illustrate now is that the experimental observation of the consequences 
of this theoretically-predicted singularity will be blurred by this unavoidable finiteness of 
the time window of any experimental observation. To see this, let us display the same 
information presented in Fig. [141( a). which plots T a (t;<f)) as a function of t for a sequence 
of volume fractions 0, in a complementary format. This is done in Fig. [TW b). which plots 
T a {t; 0) as a function of for a sequence of waiting times t. 

The main feature to notice in each of the curves corresponding to a fixed waiting time t, 
is that one can distinguish two regimes in volume fraction, namely, the low-0 (equilibrated) 
regime and the high-0 (non-equilibrated) regime, separated in a continuous fashion, and not 
as an abrupt transition, by a crossover volume fraction (f)( c >(i). Focusing, for example, on 
the results corresponding to t = 10 3 , we notice that <fyd{t = 10 3 ) ~ 0.57. In Fig. fl4T b) we 
have highlighted the crossover points <p — 0^(^) ; T a(t,<fi) = T^ q ((f>^), corresponding to each 
waiting time t considered. We observe that the resulting crossover volume fraction <fy^{t) 
first increases rather fast with t, but then slows down considerably, reaching a theoretical 
maximum crossover volume fraction, limt_ s . 0O <p^ c \t), given by (j)^ = 0.582, as indicated in 
the inset of the figure. 

The scenario illustrated by Fig. [14T b) has additional physical implications. Although it 
is impossible to witness the infinite-time implications of the theoretically-predicted singular 
dynamic arrest transition, it is important to stress that its finite-time consequences, such as 
those illustrated in this figure, can be predicted, and could be corroborated by performing 
measurements at intentionally finite, accessible waiting times. It is thus important to test if 
these predictions make sense by comparing them with available experimental or simulation 
data. Although this comparison falls out of the scope of the present paper, we can say 
that the picture that emerges from the predicted dependence of r a (t;<p) on waiting time 
and volume fraction just discussed, is fully consistent with the most relevant qualitative 
features observed in a simulation experiment consisting precisely of the equilibration of a 
hard-sphere liquid, initially prepared in a non-equilibrium state 8]. In fact, such simulation 
experiment was originally inspired by the very initial version of the theoretical scenario of 
the present non-equilibrium theory. To have an idea of the level of agreement, in the inset 
of Fig. [14T b) we have included the simulation data for the evolution of the crossover volume 
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fraction (p^ c \t) with waiting time reported in Ref 8|. In a separate communication 35] we 
shall analyze in detail other aspects of the comparison between the predictions of the NE- 
SCGLE theory and available simulation results, which indicates a general agreement and 
exhibits some well-defined limitations of the present non-equilibrium theory 



VII. CONCLUDING REMARKS 



In summary, in this work we have started the systematic exploration of the predicted NE- 
SCGLE scenario of the irreversible isochoric evolution of a soft-sphere glass-forming liquid 
whose temperature is suddenly quenched from its initial value T"W to a final value T^> = 0. 
As we explained here, the response falls in two mutually exclusive possibilities: either the 
system will reach its new equilibrium state within an equilibration time i e<? (0) that depends 
on the fixed volume fraction 0, or the system ages forever in the process of becoming a glass. 

In the first case the equilibrium a-relaxation time T^ q (4>), and the equilibration time 
t eq (4>) needed to reach thermodynamic equilibrium, are predicted to remain finite for vol- 
ume fractions smaller than a critical value 4>^ ~ 0.582, but as approaches this hard-sphere 
dynamic-arrest volume fraction, both characteristic times will diverge and will remain in- 
finite for > <p( a \ Although it is intrinsically impossible to witness the actual predicted 
divergence, the theory makes distinct predictions regarding the transient non-equilibrium 
evolution occurring within experimentally-reasonable waiting times t, which could, thus, be 
compared with realizable experiments or simulations. 

This applies even more to the predictions regarding the complementary regime, > 0* a * ) , 
in which the system, rather than ever reaching equilibrium, is predicted to age forever. As 
discussed in the previous section, under these circumstances the long-time asymptotic limit 
of S(k;t) will no longer be the expected equilibrium static structure factor S^ eq \k), but 
the non-equilibrium, but well-defined, dynamically arrested static structure factor S^(k). 
Furthermore, S(k; t) is predicted to approach S^ a \k) in a much slower fashion (a power law), 
in contrast with the exponential-like manner in which S(k;t) approaches S^ eq \k) when the 
system equilibrates. 

Putting together the two regimes just described, we have presented the scenario predicted 
to emerge for the crossover from equilibration to aging. As discussed in the previous section, 
the discontinuous and singular behavior is intrinsically unobservable in practice, due to the 
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finiteness of the time windows of experimental measurements. This forces the discontinuous 
dynamic arrest transition to appear as a blurred crossover, which may depend on the protocol 
of the experiment and of the measurements. Testing these predictions by comparing them 
with available experimental or simulation data is an issue that we shall leave for future 
studies, since the main purpose here was to provide the details of the methodologies needed to 
solve the equations that define the NE-SCGLE, and to illustrate its use with the application 
to the specific system and processes considered here. As indicated at the end of the previous 
section, we can say that the picture that emerges from the predicted dependence of r a (t; 0) 
on waiting time and volume fraction is consistent with the most relevant qualitative features 
observed in the simulation experiment of the equilibration of the hard-sphere liquid [8j. In 
a separate paper we shall establish a more direct contact with those simulation results, and 
with other simulation or experimental data. 

In the meanwhile, it will also be interesting to interrogate the NE-SCGLE theory on the 
variations of the scenario just described, when the system and conditions employed here 
are modified. For example, one may be interested in understanding how this scenario might 
change when protocol of the quench is modified. Other questions may refer to the dependence 
of this scenario on the particular class of model system and interactions (involving here only 
soft repulsions), particularly when attractive forces are incorporated. The answer to these 
questions will surely use the methods and experience developed in the presented work, and 
will be the subject of future research. 
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